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ABSTRACT 

Optimal decoding of bit interleaved coded modulation (BICM) 
MIMO-OFDM where an imperfect channel estimate is available at 
the receiver is investigated. First, by using a Bayesian approach in- 
volving the channel a posteriori density, we derive a practical de- 
coding metric for general memoryless channels that is robust to the 
presence of channel estimation errors. Then, we evaluate the out- 
age rates achieved by a decoder that uses our proposed metric. The 
performance of the proposed decoder is compared to the classical 
mismatched decoder and a theoretical decoder defined as the best de- 
coder in the presence of imperfect channel estimation. Numerical re- 
sults over Rayleigh block fading MIMO-OFDM channels show that 
the proposed decoder outperforms mismatched decoding in terms of 
bit error rate and outage capacity without introducing any additional 
complexity. 

1. INTRODUCTION 

Multiple-input multiple-output (MIMO) antennas systems is a promis- 
ing technique for high-speed, spectrally efficient and reliable wire- 
less communications. However, as higher data rates lead to wide- 
band communications, the underlying MIMO channels exhibit strong 
frequency selectivity. By using orthogonal frequency division mul- 
tiplexing (OFDM) and applying a proper cyclic prefix (CP), the fre- 
quency selective channels are transformed to an equivalent set of 
frequency-flat subchannels. These considerations motivate the com- 
bination of MIMO and OFDM, referred to as MIMO-OFDM, as 
a promising technology for the future generation of wireless sys- 
tems [1]. 

It is well known that reliable coherent data detection is not possi- 
ble unless an accurate channel state information (CSI) is available at 
the receiver. A typical scenario for wireless communication systems 
assumes the channel changing so slowly that it can be considered 
time invariant during the transmission of an entire frame. In such sit- 
uations, pilot symbol assisted modulation (PSAM) has been shown 
to be an effective solution for obtaining CSI. It consists of inserting 
known training symbols (pilot) at the beginning of the information 
frame from which the receiver estimates the channel before decod- 
ing the rest of the frame based on that channel estimate. However, 
in practical systems, due to the finite number of pilot symbols and 
noise, the receiver has only access to an imperfect (possibly very 
bad) estimate of the channel. 

In order to deal with imperfect CSI, one sub-optimal technique, 
known as mismatched decoding, consists in replacing the exact chan- 
nel by its estimate in the receiver metric. Although, this scheme is 
not adapted to the presence of channel estimation errors (CEE), it 
has been extensively adopted for performance evaluation of single an 
multi-carrier MIMO systems [2]. However, adopting a mismatched 



decoding approach arises two important questions: i) what are the 
limits of achievable information rates and how close it is to the rates 
obtained by a theoretical decoder in the presence of CEE ? and ii) 
what type of practical encoder/decoder can achieve the best perfor- 
mance under imperfect channel estimation ? The first problem has 
been addressed in [3] where mismatched decoding has been shown 
to be largely suboptimal compared to the capacity provided by a the- 
oretical decoder. Since the theoretical encoder/decoder can not be 
implemented on practical communication systems, we have recently 
addressed in [4] the second question in the case of multicarrier sys- 
tems by deriving a practical decoding metric adapted to imperfect 
channel estimation and its associated achievable rates. 

In this work, we address the above problems for more general 
wideband block fading MIMO-OFDM channels estimated by a finite 
number of training sequence. In our approach we exploit an inter- 
esting feature of PSAM, the availability of CEE statistics, in order to 
derive the a posteriori pdf of the perfect channel conditioned on its 
estimate. The latter let us to formulate a new channel likelihood as a 
function of the estimated channel and derive our improved decoding 
rule adapted to the presence of CEE. Interestingly, the present metric 
coincides with that proposed in [5] for MIMO space-time decoding. 
Furthermore, we evaluate the capacity associated to the improved 
and mismatched decoders by using the results obtained in [6]. Al- 
though in this work we will focus on wideband MIMO-OFDM sys- 
tems working in a bit interleaved coded modulation (BICM) frame- 
work, our results can be applied to general memoryless channels. 

The outline of this paper is as follows. In Section|2]we describe 
the MIMO-OFDM channel and its pilot based estimation. We also 
calculate the posterior distribution of the perfect channel conditioned 
on its estimate. This posterior distribution is used in section[3]to de- 
rive the improved decoding metric in the presence of imperfect chan- 
nel state information at the receiver (CSIR). Section|4] provides the 
application of the improved metric for iterative decoding of BICM 
MIMO-OFDM. In section|5] we evaluate the achievable outage rates 
of a receiver using the proposed metric. Section[6]illustrates via sim- 
ulations, a comparative performance study of the proposed decoder 
and section|7jconcludes the paper. 

Notational conventions are as follows. Upper and lower case 
bold symbols are used to denote matrices and vectors, respectively; 
In represents an N x N identity matrix; E x [.] refers to expectation 
with respect to x; |.| and ||.|| and Tr(.) denote matrix determinant, 
Frobenius norm and matrix trace respectively; (.) T and (.) n denote 
vector transpose and Hermitian transpose, respectively. 

2. SYSTEM MODEL AND CHANNEL ESTIMATION 

We consider a single-user MIMO-OFDM communication system 
over a memoryless frequency selective Rayleigh fading channel. The 
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Fig. 1. Block diagram of MIMO-OFDM transmission scheme. 



system consists of Mt transmit and Mr receive antennas (Mr > 
Mt), and M is the total number of subcarriers. Fig. Q] depicts the 
BICM coding scheme used at the transmitter. The binary data se- 
quence b are encoded by a non-recursive non-systematic convolu- 
tional (NRNSC) code before being interleaved by a quasi-random 
interleaver. The output bits d are multiplexed to Mt substreams and 
mapped to complex M C -QAM symbols before being modulated by 
the OFDM modulator and transmitted through Mr antennas. 

Let s be the MMt X 1 vector containing the OFDM symbols 
transmitted simultaneously over Mt antennas. The symbols are as- 
sumed to be independent identically distributed (i.i.d.) with zero 
mean and unit covariance matrix S s = E[ss H ] = Imxm t - As- 
suming an invariant channel over a frame of L symbols, the received 
vector y at a given time index I (omitted for brevity) can be written 
as 

y = Hs + 2 (1) 

where H is a MMr x MMt block diagonal channel matrix con- 
taining the frequency response of the MIMO channels and the noise 
vector z is assumed to be a zero-mean circularly symmetric com- 
plex Gaussian (ZMCSCG) random vector with covariance matrix 
S z = E(zz n ) = ct^Imxa/r- We assume that for each frame, a 
different realization of H independent of both s and z is drawn and 
remains constant during this frame. The MIMO-OFDM channel can 
be decoupled into M frequency flat MIMO channels by exploiting 
the block diagonal structure in (TJl which can be rewritten as a set of 
M equations that contains only one subcarrier each 



y k — H fc s k + z k k= 1, 



M, 



(2) 



where H = diag[HiH 2 ... H A /], y T = [yj ■ ■■ yh\, s* — [si ... s 
and z T = \z\ ... zj { ]. The architecture of ((2} constitutes the basis 
for the study in this paper. 
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Pilot Based Channel Estimation: We consider the estimation of 
channel matrix Hfe via the transmission of N training vectors s T t , 
(i — 1, N). According to (O, when pilot symbols are transmitted 
we receive 

Yt = Hfc St + Zt (3) 

where each column of the Mt X N matrix St = [s T 1 \...\s T N ] 
contains one pilot symbol and the noise Zt has the same distribu- 
tion as the noise Zk- The average energy of the training symbols is 
Pr = TvT7zTi:(StSt). The ML estimate of H fe is obtained by 



NM T \ 

minimizing ||Yt 



H fe St || 2 with respect to Hfe. We have 



til 



Yt St (StSt ) 



H fe +£: 



(4) 



where £ — ZtSt (StSt ) 1 denotes the estimation error matrix. 
When the training sequence is orthogonal (S T S^ = NP t Im t ), 
the j-th column £j of the estimation error matrix reduces to a white 



noise vector with covariance matrix £ £ = E[£ j£ "] = cr £k ^M R 
2 o- 2 

(for all j), where a £ k — -j^p— • By assuming that the channel ma- 
trix Hfe has the prior distribution Hfe ~ C/V(0, Ia/ t ® Sh.h), 
(Xif,fc = (T^Imr) and choosing an orthogonal training sequence, 
we can derive the posterior distribution of the perfect channel condi- 
tioned on the estimated channel as 



/(Hfe|Hfe" IL ) = CA/"(£ a Hml, 1m t 



SaS; 



(5) 



y £,k 



where S A = S H ,fe(E £ + En.fc)" 1 = si Mr and s = 

The availability of the estimation error distribution constitutes an 
interesting feature of pilot assisted channel estimation that we used 
to derive the posterior distribution 10. This distribution is exploited 
in the next section, in the derivation of a new metric for improving 
the detection performance under imperfect channel estimation. 

3. MAXIMUM-LIKELIHOOD DETECTION IN THE 
PRESENCE OF CHANNEL ESTIMATION ERRORS 

3.1. Mismatched ML Detection 

It is well known that under i.i.d. Gaussian noise, detecting s k is 
given by maximizing the likelihood function ^(yfelsfe, Hfe) which 
is equivalent to minimizing the Euclidean distance T>ml 

sjf^Hfc) = argmin { £>M L (sfc, y k , Hfe)}, (6) 

s fc eC M TX! 



where ©ML(sfc,yfe, Hfe 



•In W(yfe|sfe,H fe ) oc ||y fe -HfeSfe 



Since the above detection rule requires the knowledge of the per- 
fect channel matrix Hfe, one sub-optimal approach, referred to as 
mismatched detection, consists in replacing the exact channel by its 
estimate in the receiver metric as 



s k (Hfe) 



arg mm 

jtS C m t>< : 



{ \\Uk - H fc Sfe|| 2 }| 



(7) 



3.2. Improved ML Detection for imperfect CSIR 

The Bayesian framework introduced previously let us to define a new 
likelihood function W(yk\Uk, Sk) by averaging the likelihood that 
would be used if the channel were perfectly known (W (yk\Hk, s k )), 
over all realizations of the perfect channel for a given estimated 
channel state (the posterior distribution (5}). This yields 



W(y k \H k , Sfe) = E Hj |fifc [W(y k \Kk,Sk)\ Hfe] 

W(y k \H k ,s k ) /(Hfe|Hfe) dHfe. 

(8) 



/ 



Since both W(yk\Uk, Sfe) and /(Hfe|Hfe) are gaussian densities, it 
is easy to show that W{y k \H.k,s k ) = CN(n M , S M ) with [4] 



= 5Hfe Sk, 

= S 2 +5S £ ||sfc|| 2 . 



Now the ML estimate of Sfe can be formulated as 

-A1 



Sfe - (Hfc)= argmin { T> M (sk, yk, Hfe)}, 

Sfc£C M T xl 



(9) 
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with 
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Fig. 2. Block diagram of MIMO-OFDM BICM receiver. 



Notice that the channel likelihood W(yk\sk,Uk) involved in 
d 1 3 1 > is conditioned on the perfect channel Hfc of which the receiver 
has solely an estimate Hfc. In order to enhance the detection perfor- 
mance under imperfect CSIR, we propose an improved decoder that 

uses W(yk\sk, Hfc) of ([8j instead of W(yk\sk, Hfc) (mismatched 
approach), for the derivation of the coded LLRs of d 1 3 1 > . The de- 
coder accepts the LLRs of all coded bits and employs the well known 
forward-backward algorithm [8] to compute the LLRs of informa- 
tion bits, which are used for the decision. 

5. ACHIEVABLE RATES OF MEMO-OFDM UNDER 
CHANNEL ESTIMATION ERRORS 

In this section we provide the instantaneous achievable information 
rates Cm an d Cml associated to a receiver using the decoding rules 
of d 1 Ob (improved) and (mismatched), respectively. We consider 
a MIMO-OFDM channel W(y\s, H) = YllU CN{H k s k , S Zfc ), 
with S Zfc = o z k Ia/jj ■ Furthermore, we assume that the transmitter 
does not disposes of the channel estimate to perform power control. 
Thus, equal power is allocated to symbols Sk by setting S Sfc = 
PIm t . 



being the new ML decision metric under CEE. 

We note that when exact channel is available (Hfc = Hfc), the 
posterior expectation of ([8) becomes equivalent to replacing Hfc by 
Hfc in W(yk\Uk, Sfc) and consequently the the two metrics T> M 
and Dml coincides. Under near perfect CSIR, occurred either when 
oo, we have 5 — > 1, Sag — ► 0, fx M — > 
. Consequently, the improved metric T> M 



o~% — > or when N 



H,, s and S , 



behaves similarly to the classical Euclidean distance metric T>ml- 

4. ITERATIVE DECODING OF BICM MIMO-OFDM 
BASED ON IMPERFECT CSIR 

As a practical application of the general decoding metric dl It . we 
consider soft iterative decoding of BICM MIMO-OFDM under im- 
perfect CSIR. This problem has been addressed in [7] under the as- 
sumption of perfect CSIR. Here, without going into the details, we 
extend the results of [4] to MIMO-OFDM block fading channels es- 
timated by a finite number of training symbols. 

As shown in Fig. [2] the BICM receiver consists of a bunch of 
demodulator/demapper, a de-interleaver and a soft-input-soft-output 
(SISO) decoder. Let d k m be the m-th coded and interleaved bit 
(m = 1,2, ...,B = log 2 M c ) of the constellation symbol Sk at 
the the i-th transmit antenna and the k-th subcarrier. We denote by 
L{d k ,m ) the coded log-likelihood ratio (LLR) value of the bit d\ m . 
At each decoding iteration, the LLR values conditioned on the CSIR 
are given by 

L(4' m ) = log Pdem(dr = 1)|yfc ' Hfc) . (12) 
fe Pde m «' m = 0|yfc,Hfc) 

We have (see [4] and references therein) 
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where P^ cc (d l f! m ) and -Pdcc^fc" 1 ) aK extrinsic information coming 
from the SISO decoder. 



5.1. Achievable Rates Associated to the Improved Decoder 

The capacity of a general mismatched decoder assuming Gaussian 
i.i.d inputs s ~ CN(0, S s ) is given in [6] by minimizing the mutual 
information 7(5; Y) defined for MIMO-OFDM channels as [1] 

M 

I(S; Y ) = E det (Im* + TfcS^T^E,- 1 ) (14) 

fc=i 

where the minimization is done over all general channels V(y\s, T) = 
rifcli CAf(TkSk, Sfc). Since in the present work the decoding is 
performed on a single subcarrier basis, the minimization can be done 
independently for each subcarrier according to its corresponding chan- 
nel and estimate pair (Hfc, Hfc). Thus, we have to find the optimal 
Tfc G (£M R xM T an[ j covar j ance matrices Sfc = allM R so as to 
minimize I(S; Y) under the set of constraints [6] 

(ci) :Tr(E s [V(y|s,T)]) = Tr(E B [W(y|s,H)]), 

(c 2 , k ) : E Sk \Ev{V M (sk,y k \Kk)}\ < 
E Sk [Ew[©Ai(sfc,3/fc|Hfc)]], 

for k = 1, M. Applying our channel model Q}, and after some 
calculus, the above set of two constraints are obtained as 



(ci) : det ( 



T£ S T' 



(c 2 , k ) : ]|T fc + afcHfcH 2 < [|H* + a fe Hfc 



det HS S H" + S 5 



Cstfc , 



(15) 
(16) 



where 



a k = S k (SkO-j M P~\ n ,kO-^ k )[M T Sa£ M X„,kP+Xn,kO- 2 z _ k ~SkO-jj.P] ; 
Cstfc = M T Xn.h [||Hfc || 2 - ||Tfc f + i (Tr(S 2fc ) - Tr(E fc ))] [l - 
s -J$ A ».fc _ M T X n ,k] ~ , and 



(13) A„,fc 



exp 



with n = Mt — 1, and 
(-1)" 



T(-n,t) = 



71 — 1 , | 

r(0,t)-exp(-t) 



where r(0, t) is the exponential integral function. 

We consider the singular value decomposition (SVD) of = 
U fe A fc V^ with A fc = diag(Afe,i, . . . , Xk,M R )- Let D Mjfc be a di- 
agonal matrix such that -D M ,fc = U^TfcVfc, whose diagonal val- 



ues are given by the vector u = (fik.i, 



HV = V^H^Ufe and the vector h 



• • , Vk,M E ) T ■ We define 
diag(H^) T resulting of 



its diagonal and let b k = ||H fe + a fe H fc || - ci fc v,||xxfcj| - ||J 
Using the above definitions, it is easy to reformulate d 1 4b as 



2 -al(||H fc || 2 Hlh fe "^ 



M Mr 



Cm(H, H) 



E E iog 2 i 
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fc=li=l \ cr 2 (£lfc) 

subject to \\fi +ahk|| < bk- 



(17) 

According to l |17t , for each subcarrier k = 1, M, the achievable 
rates associated to the metric T>m of d 10b . can be obtained by using 
standard Lagrange method. We have 



1 Jt f 
Cm(H, H) = — J2 l°g 2 det ( Im r + 



± opt, A: - 1 opt,/c 



(18) 



ll7(l|A,|| 2 



|£V,fc|| 2 ) + a 2 Z:k , 



where ^(/x* J 
T" opl ,fc = UfcD^opi and the vector containing the diagonal ele- 

ments of the optimal solution D opt is given by 
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if b fc > 0, 



otherwise. 



(19) 



Similarly, we can compute the achievable rates associated to the 
mismatched ML decoder 10. This is given by replacing the vector 

M opl by 

_ M e {Tr(A«h fe )}~ 



(20) 



in equation <18t . 



5.2. Outage Rates Evaluation 

Given any pair of matrices (H, H), the expression d 1 8I > provides the 
instantaneous achievable rates associated to a receiver using the de- 
coding rule d 1 Ob . The outage probability associated to an outage rate 
R > is defined as 

P™\R,H) = Pr H| H(H G \ M (R,fi)\B), 

with A M (R, H) = {H : C M (H, H) < 7?}. Using this, the 
outage rate of the improved decoder for an outage probability 7 is 



G 



M ( 7 ,H)=sup{ J R>0:P^ t (ii,H)< 7 } 

R 



(21) 



and Cml'^j H) (for mismatched ML decoder), is given by replac- 
ing A^ (R, H) by A ML (i?, H) = {H : Cml(H,H) < R} in 
equation i2H . Since these outage rates still depend on the random 
channel estimate H, we will consider the expected outage rates over 
all channel estimates as 



C7XT(7) = E ft [07,H) 



(22) 



The achievable rates i22l are upper bounded by the outage rates pro- 
vided by a theoretical decoder (i.e. the best decoder in the presence 
of CEE). In our case, this capacity is given by 



C; ut ( 7 )=EH[cr(7,H)] 



where the outage rates C° ut (7, H) are computed using 



g 

M 



C s (H,H) = ] if>g2 det (Im r + 



P H fe H^ 



(23) 



(24) 



and Hfe are random channels drawn from the posterior distribution 
of©. 

6. NUMERICAL RESULTS 

Next, the performance of the improved decoder is measured in terms 
of bit error rate (BER) and achievable outage rates. The binary infor- 
mation data are encoded by a rate 1/2 NRNSC code with constraint 
length 3 defined in octal form by (5,7). Throughout the simulations, 
each frame is assumed to consists of one OFDM symbol with 50 
subcarriers belonging to a 16-QAM constellation with Gray label- 
ing. The interleaver is a pseudo-random one operating over the en- 
tire frame with size M ■ Mt ■ log 2 (M c ) bits. For each transmitted 
frame, a different realization of a Rayleigh distributed channel has 
been drawn and remains constant during the whole frame. Besides, 
it is assumed that the average pilot symbol energy is equal to the 
average data symbol energy. Moreover, the number of decoding it- 
erations are set to 4. 

Fig. [3]depicts the BER performance over a 2 x 2 MIMO-OFDM 
channel estimated by N £ {2, 4, 8} pilot symbols per frame. As 
observed, the increase in the required Eb/No caused by CEE is an 
important effect of imperfect CSIR in the case of mismatched ML 
decoding. The figure shows that the SNR to obtain a BER of 10 
with N — 2 pilots is reduced by about 1.3 dB if the improved de- 
coder is used instead of the mismatched decoder. We also notice that 
the performance loss of the mismatched receiver with respect to the 
derived receiver becomes insignificant for N > 8. This can be ex- 
plained from the expression of the metric i ll It , where we note that by 
increasing the number of pilot symbols, this expression tends to the 
classical Euclidean distance metric. However, the proposed decoder 
outperforms the mismatched decoder especially when few numbers 
of pilot symbols are dedicated for channel estimation. 

Fig. E] shows the expected outage rates (in bits per channel 
use) corresponding to the transmission of one OFDM symbol with 
M = 16 subcarriers and Mt = Mr = 2 antennas, achieved by 
adopting mismatched ML decoding and the improved decoder (ex- 
pression l!22t). For comparison, we also display the upper bounds 
on these achievable outage rates (expression l |23t ) and the ergodic 
capacity. N — 2 pilot symbols are sent per frame for CSIR acqui- 
sition and the outage probability has been fixed to 7 = 0.01. At a 
mean outage rate of 8 bits, we note that the achievable rate of the 
mismatched decoder is about 5 dB of SNR far from the capacity 
provided by the theoretical decoder. As observed, by using the im- 
proved decoder, higher rates are obtained for any SNR values and 
the aforementioned SNR gap is reduced by about 1.8 dB. 

Similarly, Fig. [5] compares the outage rates in the case of a 
4x4 MIMO-OFDM channel estimated by N = 4 training sym- 
bols. Again, it can be observed that the proposed decoder achieves 
higher rates than the mismatched decoder. However, we note that 
the increase in the SNR (at a given mean outage rate) induced by 
using the mismatched decoder rather than the improved decoder, is 



less than that obtained for a 2 x 2 MIMO-OFDM channel (see Fig. 
0. This can be explained by noting that when Mr = 4, we must 
send N > Mt = 4 pilot symbols per frame which yields a more 
accurate estimate of the channel and consequently bring closer the 
performance of the two decoders. This observation is in consistence 
with those presented in [2] where it is reported that the performance 
degradation due to imperfect channel estimation can be reduced by 
increasing the number of antennas. 

7. CONCLUSION 

A Bayesian approach for the design of a decoding method that is 
robust to channel estimation errors was presented. The robustness 
of our method comes from the averaging of the decoding rule, that 
would be used if the channel were perfectly known, over all chan- 
nel estimation errors. We also derived the expression of the achiev- 
able rates associated to our proposed decoder and compared it to 
that provided by the classical mismatched decoding approach. As 
a practical application, the proposed decoder was exploited for it- 
erative BICM decoding of MIMO-OFDM under imperfect channel 
knowledge. Simulation results showed that, without introducing any 
additional complexity, the proposed decoder outperforms the clas- 
sical mismatched approach in terms of BER and achievable outage 
rates for short training sequences. Although our proposed method 
outperforms mismatched decoding, the derivation of a practical de- 
coder achieving the maximum outage rate under imperfect channel 
estimation is still an open problem. 
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Fig. 3. BER performance of the proposed and mismatched decoders 
in the case of a 2 x 2 MIMO-OFDM Rayleigh fading channel with 
M = 50 subcarriers for training sequence length N 6 {2, 4, 8}. 
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Fig. 4. Expected outage rates of a 2 x 2 MIMO-OFDM system with 
M — 16 subcarriers versus SNR for N — 2 pilots. 
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Fig. 5. Expected outage rates of a 4 x 4 MIMO-OFDM system with 
M — 16 subcarriers versus SNR for TV = 4 pilots. 



